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Abstract — For the K user, single input single output (SISO), 
frequency selective interference channel, a new low complexity 
transmit beamforming design that improves the achievable sum 
rate is presented. Jointly employing the interference alignment 
(IA) scheme presented by Cadambe and Jafar in |1| and linear 
minimum mean square error (MMSE) decoding at the trans- 
mitters and receivers, respectively, the new IA precoding design 
improves the average sum rate while preserving the achievable 
degrees of freedom of the Cadambe and Jafar scheme, K/2. 

Index Terms — Interference alignment, interference channel, 
degrees of freedom, coordinated interference mitigation, SISO, 
sum rate. 



I. Introduction 

Cadambe and Jafar (CJ) have shown in [ 1 1 that for the K 
user, single input single output (SISO), frequency selective 
interference channel, K/2 degrees of freedom is achievable 
by using interference alignment (IA) at the transmitters and 
zero forcing (ZF) decoding at the receivers. Recent work by 
Shen, Host-Madsen, and Vidal (SHV) has shown that the joint 
IA and ZF strategy can be further improved to increase the data 
rate performance of the overall system while maintaining the 
achievable degrees of freedom in the following two ways. One 
in which new precoding subspaces are used to generate the 
precoding vectors and another which optimizes the precoding 
vectors within the subspaces defined in the original CJ scheme. 
While the optimization of the second method is presented in 
0, an explicit solution to the first is not. 

Then in [3|, Kim and Torlak (KT) proved that the sum 
rate of the K = 3 SISO interference alignment channel in 
which receivers employ linear minimum mean square error 
(MMSE) decoders is a concave function with respect to w, 
(fl2] i in Sec. II, a vector with which the precoding subspaces 
are defined. Moreover, they prove that the concave sum rate 
expression can be formulated into a constrained maximization 
problem for which a global solution exists thereby improving 
the conditioning of the precoding subspaces and maximizing 
the sum rate with respect to w. 

However, calculation of the optimal solution requires use 
of a first order optimization algorithm, i.e. gradient descent. 
Therefore, as a trade off in data rate performance for reduced 
computational complexity, KT also derive a lower bound to the 
sum rate expression that when maximized yields a solution that 
is far less complex to calculate but still achieves significant 



increases in sum rate while also maintaining the achievable 
degrees of freedom of the CJ scheme. 

In this paper, a new lower bound to the sum rate expression 
introduced in is derived that when maximized produces a 
simple closed-form solution that greatly reduces computational 
costs compared to either of the solutions proposed by KT while 
still providing considerable gains in data rate performance. 
Moreover, the new closed-form solution also maintains the 
achievable degrees of freedom proposed by CJ. 

Notations: we use upper case letters to denote matrices, 
lower case letters to denote scalars, and boldface letters to 
denote vectors. (•)* and {-y refers to the transpose and 
conjugate transpose of (•), respectively. 

II. System Model 

Consider the K = 3 user interference channel where user 
1 is allocated n + 1 streams of data while users 2 and 3 are 
both allocated n; and all three users transmit over a 2n + 1 
symbol extension of the channel by means of beamforming. 
Here, n € N. The channel output at the kt\\ receiver is then 
defined as 

y k := H kl ViXi + H k2 V 2 x 2 + H k3 V s x s + z k 

where k — {1,2,3} is the user index. y k is an N by 1 vector 
representing the output signal of the kth receiver where N = 
2n + 1. H k j is a diagonal N by N matrix of channel fading 
coefficients from transmitter j to receiver k. 

The (n + 1) by 1 vector X\ represents the input signal 
of the first transmitter and x k for k = {2, 3} is an n by 
1 vector representing the input signal of the kth transmitter. 
Vector Xi ~ CN(0,pI n+ i), and x 2 , x$ are both CN(0,pI n ). 

Precoding matrix V\ is N by (n + 1) and matrices V 2 and 
V3 are both N by n. 

Lastly, z k is an N by 1 vector representing the additive 
white Gaussian noise (AWGN), i.e. z k ~ CN(0, I N ). 

In order to guarantee the alignment of interference at each 
of the three receivers, CJ define the matrices of beamforming 
vectors as 
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and the N by 1 vector w for the CJ scheme is a vector of all 
ones. 

Once aligned, the interference at receivers 1, 2, and 3 are 
constrained respectively to the following. 



H12V2 = H13V3 
H23V3 -< H 2 \Vi 
H32V2 -< H31V1 



(8) 
(9) 
(10) 



where A -< B denotes the set of column vectors of matrix A 
is a subset of the column vectors of matrix B. More precisely, 
the set of columns in H23V3 are equal to the last set of n 
columns in matrix H21V1. Likewise, the set of columns in 
H32V2 are equal to the first set of n columns in matrix H31V1. 
For example, if H32V2 — H31C, then 

H31V1 =H 3 i[C T n w] = [ H32V2 H 31 T n w ] . 

In 0, the three beamforming matrices are redefined as 

V k (w):=W(w)T k (11) 

where the N by N diagonal matrix, W (w), is defined as a 
function of the TV by 1 vector 



w = [ W\ W2 



W N 



(12) 



such that (W (w)) u := Wi and Wi € R + for all i. Note that by 
setting w equal to 1, the new expression for the beamforming 
matrices in (fTTb is equivalent to the set of precoders for the 
CJ scheme. But by factoring out the diagonal elements of W 
from Vk and defining them as variable, a new w can be found 
that improves the sum rate performance. 
The N by n + 1 matrix Ti is defined as 

r x := [ 1 ti ■■■ t n ] 

and the N by n matrices 1?2 and are respectively defined 

as 

T 2 := (I^r 1 #31 [1 ti ••• t„_i ] 
^3 '■= {H23) 1 H 2 i [ t\ t 2 ■■■ t n 

where t m — diag(T m ) for m = {l,...,n}. The function 
diag(-), like that within Matlab, creates a column vector 
comprised of the diagonal elements of its matrix input. 

III. Achievable Sum Rate and Lower Bound 
Approximation 

In this section, we introduce the sum rate expression derived 
in Q with which we formulate a new lower bound. 



Introduced in yj, the sum rate for the K = 3 SISO 
interference channel in which receivers employ linear MMSE 
decoders is defined as 
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where the individual rate for user k is defined as 
det (i + p J2j IliAVVjV^V'll, 



fk(w) := log 



dot (/ + p E j¥fc H^WTjT^Wml 



(14) 



for j = {1,2,3}. Note that the dependence on vector w in 
(fl~4T > is implied to be within diagonal matrix W as in (fTTb 
though not explicitly expressed here nor for the remainder of 
this paper as for the sake of brevity. 

In Lemma 1 of 0, equation dT3b is expressed equivalently 
in the following compact form 



- log dot (/ + P G{ d WG kd ^ (15) 

where W := W*W and w is the vector formed of the N 
diagonal elements of W, i.e. Wi :— \wi\ 2 . Furthermore, 



and 



Gin 

Gid 

G2n 
G2d 
G3 n 
G3d 



P 2 ■= 
P3 ■ = 



[ HxxYi Gid 
V2H 12 T 2 

[ H 2 2^2 
H21T1P2 
[ H33T 3 

H 3 ir 1 P 3 



G2d 

G 3 d 



(16) 
(17) 
(18) 
(19) 
(20) 
(21) 



1 

Onxl 

Olxn 



Olxn 

V2I n 

Onxl 
1 



It is based upon this expression of the sum rate which we 
derive our new lower bound. 



B. Achievable Sum Rate Lower Bound 

In deriving the lower bound to ( TTBI l. we begin by taking its 
limit as p goes to infinity as shown in (1221 at the top of the 
next page. 

Next, we derive an upper bound to the magnitude of the 
final negative term in d22l . denoted as /„{, (w), that in turn 
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lower bounds the high SNR approximation of the sum rate. 
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where d23l applies Hadamard's inequality [4] and d24l > follows 
from x > logx for x > 0. Vector g fedi is the ith column of 
Gkd- 

Next, we express the sum of traces in (l25l l as a trace of 
sums, 



f ub (w) = ^Tr 



^E G w Gfcd 



and substitute the interference alignment constraints (HJ - dTol > 
into ([T7]i, ([T9]l, and (ED to arrive at 
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where 7 fci is the ith row vector of matrix T k and 

|2 



£l(*i*)< 



for fc = {1, 2, 3}, j = {1, 2, 3}, and i = {1, 2, TV}. 

We then upper bound (|26| | by replacing the coefficients c te i 
with the single coefficient c such that c := maxfcn} for all 
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and i. The final expression of the upper bound f u b (w) is 
en 



fub (*) = ]y E ^ 



^E r * r i 



Thus, the lower bound to the sum rate ( fT3T > is defined as 
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Note that the slope of the first term on the righthand side of 
the above states that the lower bound, as a function of logp, 
increases linearly with a slope of (3n + l)/(2n + 1) as is 
the case with (T5[ when the limit of p approaches infinity, 
thus maintaining the achievable degrees of freedom of the CJ 
scheme. 

IV. Suboptimal Improved Subspace Design 

Based upon the lower bound to the sum rate d27| >, we present 
the main result of this paper summarized in the following 
theorem. 

Theorem 1: The maximization of the sum rate lower bound 
defined in d27l i for the K = 3 SISO frequency selective 
interference channel is defined as 



C'w = max fi b (w) 



subject to 



E Tr [ 



nwr k 



= 3 TV, Wi > for all i 



k=l 



for which the closed-form solution is defined as 

3 E^-tL 



W; 



(28) 



(29) 



for all fc = {1, 2, 3} and i = {1, 2, TV}. 

While an optimal solution to the maximization of the 
complete sum rate < fT~5T > is proven to exist in O, it nonetheless 
must be computed using a constrained numerical optimization 



routine. Another suboptimal solution to the constrained sum 
rate maximization problem is also proposed in [3 1 that while 
less complex to solve than the optimal solution, is still non- 
trivial. However, the closed form solution to the maximization 
problem above while suboptimal to those proposed in [3 1 is far 
less complex to compute and still achieves significant gains in 
data rate performance as will be demonstrated in the upcoming 
section. 

To derive the solution ( 1291 ). we define the Lagrangian 
and gradient for the above constrained maximization problem 
respectively as 

L (w, A) := f lb (w) - X (j2 Tr [ r I^ r fe] - 3N ) 

and 

V*£ (w, A) := Idiag (W" 1 ) diag f £ r fc rA 



we set the gradient above 



Next, using the constraint 
to zero and multiply it by w 



#' V A L (w, A) = 3 - y— + A J 3N = 0. 

Solving for A, we get A = (l — c)/N. Then setting the gradient 
to zero, substituting in the above solution for A, and solving 
for w, we arrive at 



V. Simulation Results 

In this section, we provide numerical results to compare the 
data rate performance of the proposed design, those presented 
in 0, and the original CJ schemeQ Data is presented for the 
following algorithms abbreviated as 

• CJ - the original CJ scheme 

• KT op - the optimal design proposed by KT in Q 

• KT sop i - the first suboptimal design proposed by KT in 

ED 

• KT sop2 - the suboptimal design proposed herein. 

Also provided as a reference is the line Dm log 2 (snr) where 
D N = {§, if} for N = {3, 11}, respectively. 

It was shown in Q that while significant gains in sum 
rate were achievable using the proposed algorithms therein, 
even further gains were achievable when combined with 
SHV orthonormalization [2|. Therefore, the simulation results 
presented for the schemes listed above, with the exception 
of CJ, have been obtained in combination with SHV or- 
thonormalization. Summaries of KT op , KT sop i, and SHV 
orthonormalization are available in the Appendix. 

As shown in Fig. Q] for N — 3, the optimal KT op design 
achieves a gain in sum rate of 2.39 bits/s/Hz over the original 
CJ scheme in the high SNR regime at 50dB. Just beneath that, 
the first suboptimal design, KT sop i, is only 0.22 bits/s/Hz less 

'Channel coefficients are drawn i.i.d. from a complex Gaussian distribution. 
However, the channel coefficients were further bounded to be within a non- 
zero minimum value and a finite maximum value in order to avoid degenerate 
channel conditions. 




Fig. 1. Average sum rate comparison under various schemes for N = 3. 
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Fig. 2. Average sum rate comparison under various schemes for N = 11. 



than the optimal design; and the proposed suboptimal closed- 
form solution, KT sop 2, is only 0.47 bits/s/Hz less than the 
optimal design. 

In the mid SNR regime at lOdB, the difference between 
KT op and KT sop 2 is 0.23 bits/s/Hz, approximately half the 
difference of that at 50dB. Moreover, a 0.73 bits/s/Hz im- 
provement in sum rate over the CJ scheme is achievable with 
use of the simpler KT sop2 design. 

Similar improvements are achievable for larger symbol 
extensions of the channel. As shown in Fig. [2] for N = 11 at 
an SNR of 50dB, the difference between KT op and KT sop2 
is 1 .04 bits/s/Hz while the difference between KT sop 2 and the 
CJ scheme is 17.56 bits/s/Hz. At a mid SNR of lOdB, the 
difference between KT op and KT sop 2 is 0.18 bits/s/Hz while 
the difference between KT sop 2 and the CJ scheme is 3.02 
bits/s/Hz. 



VI. Conclusion 

In this paper, we have presented the design of a new low 
complexity algorithm that improves the sum rate performance 
of the Cadambe and Jafar [ 1 1 K — 3 user interference 
alignment scheme while preserving its achievable degrees of 
freedom. The proposed design is based upon the subspace opti- 
mization algorithm introduced by Kim and Torlak in 1 3 ] which 
maximizes the network sum rate with respect to w. Unlike the 
optimal solution in [3], the proposed design does not require 
use of a numerical optimization solver. Rather, the proposed 
design has a closed-form solution that, while suboptimal, still 
achieves significant gains in data rate performance at far less of 
a computational cost. The proposed algorithm can be extended 
to a greater number of users and to other variations of the 
original CJ scheme Q (6). 

Appendix A 

Summary of KT Improved Subspace Optimization 

In Q, the sum rate dBJl for the K = 3 user SISO 
Gaussian constant interference channel is proven to be a 
concave function with respect to w that can be expressed in 
the following optimization problem as 

C = min minmax [/ (w) - A (q l w - 3JV)1 (30) 

q:qi<c, A>0 U>i~0 lJ V ' V J1 

where the elements of the N by 1 vectors q and c are defined 
respectively as 



Ci - Xi/X 



E 



and j ki is the ith row of IV Thus there exists a maximum C 
for some w = d* which satisfies 

— ^(Aki-Bu) = A*, if m > 

C * k 

- J2( A ki -Bki) < A*, if wi = 

A k 

for some A* > 0. Terms Aui and Bi-i are defined as 
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where g kni and g kdi are the ith row vectors of matrices Gkn 
and Gkd, respectively. 

Appendix B 

Summary of KT Suboptimal Subspace Conditioning 

In 0, a suboptimal maximization problem to that of (l30l 
is defined as 

Cm = max ./jf, (w) 



subject to 



£Tr[rt 



= ZN, m>0 for all i 



(31) 



where a lower bound to the sum rate (15[ is defined asymp- 
totically in SNR as 

3n + 1 3 
fib (w) := — logp + — log det W 



The solution to the above maximization problem is 

Wi = ( Xi b E flfci + E bkiClki ) 

V k k / 

where 
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for j = {1, 2, 3} and must satisfy 

J2k bkiClki 



A/6 > 



J2k a ki 



for all i. 



Appendix C 
Summary of SHV Orthonormalization 

In 0, for the K = 3 SISO interference alignment channel, 
Shen, Host-Madsen, and Vidal show that the high SNR offset 
of the sum rate, eq. (5) in Q, can be maximized under 
constraints (8) - (20) in (2| by precoding the CJ transmit 
beamforming vectors V2 and V3 by the n by n matrices F 
and E, respectively, such that 



r 2 



V 2 F 
V 3 E 



k=l 



and matrix V\ is unaltered, defined as in the original CJ 
scheme. 

They show that the solution to F and E orthogonalizes the 
n columns of V2 and V3, respectively, while simultaneously 
allocating N/n units of power to each column. SHV also 
show that by doing so, the solutions uphold the interference 
alignment constraints (b) in [2], maintaining the achievable 
degrees of freedom proposed by CJ. 
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